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1.4 GHz and 60 /im luminosities 
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ABSTRACT 

We determine and examine the star-formation rates of 50 well known, local Blue Compact 
Dwarf (BCD) galaxies based on their 1.4 GHz and 60 /im luminosities. We find that in cases for 
which both radio and far-infrared luminosities are available, the resulting star-formation rates 
agree extremely well with one another. We determine that the star-formation rates of the BCD 
galaxies in our sample span nearly five orders of magnitude, from approximately a few times 10"'^ 
to several times 10^ M0 yr"-'^, with a median SFR of about 0.3 M© yr^-'^. We discuss trends of 
metallicity (primarily oxygen abundance) with star-formation rate, and explore the connections 
between SFR and galaxy mass estimates. 

Subject headings: galaxies: dwarf — galaxies: starburst — galaxies: evolution 



1. Introduction 

Blue Compact Dwarf (BCD) galaxies pro- 
vide evidence that isolated, small, low-metallicity 
galaxies may experience vigorous star formation 
in the present epoch, challenging our current un- 
derstanding of galaxy formation. The very na- 
ture of BCDs remains puzzling in spite of thirty 
years of research since they were first recognized 
as an interesting class of galaxy (Sargent & Searle 
1970). For instance, it is still debated whether 
some "primeval" BCDs are perhaps in the process 
of undergoing their very first starburst, or whether 
BCDs undergo several episodic starbursts sepa- 
rated by long quiescent intervals, what triggers the 
star formation in BCDs, and whether starburst- 
induced supernova-driven winds can cause the 
transformation of BCDs into one of the other 
known types of dwarf galaxy (Kunth & Ostlin 
2000, for a recent review). While star- formation 
rates (SFRs) are central to discussions of the evo- 
lution of BCDs, there are surprisingly few quanti- 
tative SFRs for BCDs available in the literature. 
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The earliest and most quoted quantitative 
SFRs are based on direct measures of the ultra- 
violet luminosity of massive stars in formation. 
Fanelli et al. (1988) modeled the International Ul- 
traviolet Explorer spectra of seven BCDs. The 
implied SFRs for a Salpeter initial mass function 
(IMF) with index 2.35 and mass cutoffs at 0.1 and 
100 M0 (which we adopt for all measurements 
given in this paper) range from about 2.5 x 10"'^ 
to 25 Mq yr-i. 

Subsequent efforts preferentially used the re- 
processed ultraviolet radiation of massive stars as 
measures for BCD SFRs, derived from either dust, 
or ionized gas emission. Sage et al. (1992) inves- 
tigated a sample of 15 BCDs, and derived SFRs 
based on far-infrared (FIR) and Ha luminosity. 
Their values range from about 0.1 to 30Moyr~^. 
Popescu et al. (1999) studied the SFRs of a sam- 
ple of 48 compact star-forming dwarf galaxies in 
voids. Using spectroscopically derived H^ lumi- 
nosities, they determined values that range from 
0.02 to 2.2MQyr-\ with a mean SFR of about 
0.5MQyr-^ Most recently, Izotova et al. (2000) 
derived the SFRs of 27 BCDs from the Second 
Byurakan Survey (SBS). Their SFRs are based on 
FIR (60 /im) and radio continuum (1.4 GHz) lumi- 
nosities. The SFRs derived from FIR luminosity 
span the range from about 0.9 to 80 Mg yr'^ Ra- 
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dio luminosities were converted to SFRs assum- 
ing that the continuum is due to thermal, free- 
free emission. These SFRs range from approxi- 
mately 2.5 to 27OM0yr~^, with a mean of about 
ISMoyr-i. 

The values available in the literature suggest 
that the SFRs of BCDs might vary over five orders 
of magnitude, from 10^"^ to lO^MQyr"^. How- 
ever, the mean SFRs given in the existing investi- 
gations are noted to also vary over two orders of 
magnitude, from O.lMQyr"^ to lOM0yr~^. It is 
therefore not clear whether and to what extent the 
wide range of BCD SFRs results from the use of 
different techniques, assumptions, or sample selec- 
tions, employed in the available studies. Clearly, 
a systematic exploration of the parameter space 
of BCD SFRs is both warranted and timely. Ac- 
curate SFRs are needed before any correlations of 
the SFR with other physical properties can be dis- 
cussed. 

This paper is meant to be a first step toward 
establishing a baseline of systematic, quantitative 
SFRs for BCDs. It was motivated by our desire 
to measure the SFRs of BCDs using a uniform 
approach and by employing techniques that are 
prone to minimal systematic errors. Specifically, 
as star-formation may occur throughout BCDs 
and not just in their very active centers, we ex- 
cluded spectroscopic methods. These require cor- 
rections for the (unknown) amount of flux missed 
outside of the aperture. We also wanted to use 
a measurement that is free from the potentially 
large and uncertain corrections imposed by inter- 
nal and foregound dust extinctions. This imme- 
diately also ruled out ultraviolet continuum and 
optical emission line indicators of the SFR. We 
therefore decided to examine the SFRs of BCDs 
based on their radio continuum and FIR luminosi- 
ties. Fortunately, public databases of FIR and 
radio flux densities that include BCDs are read- 
ily at hand. For completeness, and to investigate 
the extent of and role played by obscuration in 
our sample, a subsequent study to be presented in 
Paper II will explore the SFRs derived from very 
wide-slit Ha spectroscopy. 

2. Sample selection 

The sample under investigation consists of 50 
of the most well known and well studied BCDs. 



It is in this regard that our sample consists of 
"typical" BCDs, and that the resulting SFRs are 
called "typical" . By no means does this investiga- 
tion present an unbiased or a complete sample of 
BCD SFRs. The present sample comes from those 
BCDs that have been well studied because they 
were used to derive the primordial helium abun- 
dance, and have been taken from a recent compi- 
lation of 54 objects (Izotov & Thuan 1999, here- 
after IT). For all of these galaxies, accurate red- 
shifts, needed to calculate luminosities, are avail- 
able. The BCDs in this sample also have care- 
ful determinations of the oxygen abundance of the 
ionized gas. Since oxygen is produced by massive 
stars, and since massive stars are used as a quanti- 
tative measure of the SFR, a particular emphasis 
of the present paper is to investigate any correla- 
tions between these two quantities. 

Observational information for the galaxies in 
this study is presented in Table 1. This con- 
tains positions ( J2000) , 60 /xm flux densities 
from IRAS (Infrared Astronomical Satellite) cat- 
alogues, 1.4 GHz flux densities from the NVSS 
(NRAO VLA Sky Survey) and FIRST (Faint 
Images of the Radio Sky at Twenty centime- 
ters) catalogues and images, and heliocentric 
velocities. Much of this information was ob- 
tained from NED^. Four of the galaxies in Ta- 
ble 1, Mrk 71, SBS 1319-h579, SBS 1331+493 and 
SBS 1533+574, were split into multiple compo- 
nents for their metallicity measurements (Izotov & 
Thuan 1999, their Table 1, which has 54 entries). 
Since the 60 /^m (IRAS) and 1.4 GHz (NVSS) flux 
densities do not have fine enough resolution to dis- 
tinguish between these components, these galaxies 
are treated as single objects in the current study. 
The metallicitics (shown in Table 2) for those ob- 
jects are the average of the values measured for 
the individual components (IT), while the errors 
reflect the ranges in values. 

2.1. Distances and flux densities 

The majority of the 50 BCDs in the sam- 
ple had their distances estimated based on he- 
liocentric velocities taken from NED. The veloc- 



^The NASA/IPAC Extragalactic Database (NED) is oper- 
ated by the Jet Propulsion Laboratory, California Institute 
of Technology, under contract with the National Aeronau- 
tics and Space Administration. 
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itics given in NED were derived from emission 
line or HI measurements made predominantly by 
Izotov, Thuan and collaborators (Thuan et al. 
1999a; Izotov & Thuan 1998; Izotov et al. 1997, 
1994), or taken from the RC3 (de Vaucouleurs 
et al. 1991). In cases where velocities for several 
components were listed, we assumed the velocity 
that is associated with the main body or main 
component of the galaxy. Distances for all but 
three of our sample were estimated from veloci- 
ties by using a Virgocentric flow model (Kraan- 
Korteweg 1986). We adopted a Hubble constant 
of 72 ± 8kms-iMpc-i (Freedman et al. 2001) 
and a Virgo distance of 17.0 ± 0.3 Mpc (Tonry 
et al. 2001). The three exceptions, where dis- 
tances were estimated from analysis of the I-band 
brightness of tip of the red giant branch (TRGB) 
stars in color-magnitude diagrams (CMDs) of stel- 
lar populations, are UGC 04483 (Dolphin ct al. 
2001), VII Zw 403 (Schulte-Ladbeck et al. 1999) 
and Mrk 0209 (I Zw 36) (Schulte-Ladbeck et al. 
2001a). We present histograms showing the dis- 
tributions of distances and heliocentric velocities 
for our BCD sample in Figure 1. 

The ratios of the distances from the Virgo- 
centric flow model to the TRGB distances for 
UGC 04483, VII Zw 403 and Mrk 0209 are 1.69, 
0.23, and <1.2, respectively. This indicates a con- 
siderable misjudgmcnt of the distances based on 
velocities and the Virgocentric flow model is pos- 
sible. We believe that distance errors can poten- 
tially amount to up to a factor of a few, but prob- 
ably do not exceed an order of magnitude. It is 
this difficulty in assessing distance errors that led 
us to omit errors on the SFRs derived below. 

2.1.1. IRAS 

The Infrared Astronomical Satellite has proved 
an immensely valuable resource in exploring the 
infrared properties of galaxies since its launch in 
1983. IRAS products include 60 /im flux densi- 
ties for galaxies from the Point Source Catalog 
(brighter than about 0.5 Jy) and Faint Source Cat- 
alog (typically brighter than 0.2 Jy). These fluxes 
were obtained for the BCD sample, where avail- 
able, from NED and from the IRSA^ ISSA image 

®The NASA/IPAC Infrared Science Archive is operated by 
the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and 



server and catalog overlay service. The resolution 
of the IRAS images at 60 /xm is about 1' with typ- 
ical positional accuracies better than 20". Of the 
50 BCDs in our sample 18 arc detected at 60 /im by 
IRAS. We have not made direct measurements for 
the flux density upper limits of the non-detections, 
but for most of these the 0.2 Jy limit of the Faint 
Source Catalog should be indicative. 

2.1.2. NVSS 

The NRAO VLA Sky Survey (Condon ct al. 
1998) is a 1.4 GHz radio continuum survey cover- 
ing the sky north of —40° declination, made us- 
ing the Very Large Array (VLA) operated by the 
National Radio Astronomy Observatory (NRAO). 
The observations were made using the compact 
D and DnC configurations of the VLA. The re- 
sulting images have an rms noise level of about 
0.5 mJy beam" ^ and a FWHM resolution of 45". 
Positional accuracies (from the source measure- 
ments using a gaussian fitting algorithm) range 
from about 1" for sources stronger than 15mJy to 
7" at the survey limit. The catalog lists sources de- 
tected above a peak brightness of 2mJybeam~^. 
For 13 of the galaxies in our BCD sample an NVSS 
catalog source is close enough (< 30") to be a 
probable identification, and the fiux measurement 
from the catalog was taken. For another 10 of 
the sample galaxies, where the NVSS FITS image 
showed a local peak in intensity at the position of 
the galaxy, but no catalog source was present, we 
performed our own source measurements. We used 
the MIRIAD (Multichannel Image Reconstruction, 
Image Analysis and Display) task sfind (Hopkins 
et al. 2002) to make these measurements, and the 
sources so detected are typically at the 4 — So- 
level. For non-detections, we measured the local 
rms noise level in the vicinity of the BCD galaxy 
and adopted 5 times this value as the upper limit 
on the NVSS flux density. Uncertainties in the 
flux density measurements of the NVSS catalogue 
have been described in detail (Condon et al. 1998, 
see also their Figure 31). In practice, these uncer- 
tainties rise from w 5% at flux densities of lOmJy 
to ~ 30% close to the survey limit. These un- 
certainties translate to uncertainties in estimated 
star- formation rates (see Section 3). Digitized Sky 

Space Administration. See http://irsa.ipa<;. caltech.edu/ 
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Survey^ images have been overlayed with 1.4 GHz 
contours from NVSS for each of the 23 detected 
BCDs (Figure 2). 

2.1.3. FIRST 

FIRST, Faint Images of the Radio Sky at 
Twenty centimeters (Becker et al. 1995), is a 
1.4 GHz radio continuum survey currently cover- 
ing about 8500 square degrees of the North Galac- 
tic Cap. The observations were taken using the 
VLA in B configuration. The images have a typ- 
ical rms noise level of 0.15 mJy with a FWHM 
resolution of 5". Source positional accuracies vary 
from 0.5" at the 3mJy level to about 1" for the 
faintest sources. The source detection threshold of 
the survey is 1 mJy. Because of the much higher 
resolution of FIRST compared to NVSS, some 
unresolved NVSS objects are resolved by FIRST 
into more than one component (UM 311 is one 
example). In this case we have taken the sum of 
the FIRST flux densities for the individual com- 
ponents as the FIRST measurement for the whole 
source. From the FIRST catalog, 6 galaxies were 
matched with our BCD sample, and we obtained 
upper limits from the FITS images for another 
14. The flux density upper limits arc again taken 
to be 5 times the local rms value. The remaining 
30 BCDs lie outside the area of the FIRST sur- 
vey. 1.4 GHz contours from FIRST are shown on 
the images in Figure 2 for each of the 6 detected 
BCDs. 

An important advantage from the higher reso- 
lution of FIRST images is being able to identify 
the positional origin of the radio emission with 
greater accuracy than allowed by NVSS. The dif- 
ferent resolutions are similarly an important effect 
to account for if calculating SFR per unit area, 
and also allow for distinguishing between the ob- 
ject of interest and confusing background sources. 
In the case of Mrk 0209, for example, the FIRST 
emission appears to originate slightly to the SE of 
the region of brightest optical emission. This may 



^The Digitized Sky Surveys were produced at the Space 
Telescope Science Institute under U.S. Government grant 
NAG W-2166. The images of these surveys are based on 
photographic data obtained using the Oschin Schmidt Tele- 
scope on Palomar Mountain and the UK Schmidt Tele- 
scope. The plates were processed into the present com- 
pressed digital form with the permission of these institu- 
tions. 



initially suggest a background source for the radio 
emission, but the strong consistency between the 
radio-derived SFR and numerous other optically- 
derived values for this galaxy, (Schulte-Ladbeck et 
al. 2001a, see also discussion in Section 5.1), ar- 
gue that in fact this radio source is likely to be 
a star-forming knot in the outskirts of Mrk 0209. 
Similarly, CG 0798 shows a double-peak in the 
NVSS emission, and a FIRST source about an ar- 
cminute to the south of the galaxy accounts for 
some of this emission. Lacking more detailed ob- 
servations, it is certainly possible that this may 
be a coincidental background source and the ra- 
dio flux-density and the derived SFR for CG 0798 
should strictly be treated as upper limits. Given 
the good agreement between the radio and FIR 
SFRs for this galaxy, however, we have chosen to 
treat it as a detection for the current analysis. 

3. Star-formation rates 

To avoid the effects of both spectroscopic flux 
losses and uncertain obscuration corrections, we 
use 1.4 GHz and 60y^m luminosities as estimators 
of the star-formation rate for galaxies in our BCD 
sample. Here we adopt the calibrations 



SFR, 



■60jum 



= 5.5 X 



5.1 X 1023(WHz-i) 



(1) 



(Cram et al. 1998; Condon 1992) for the SFR de- 
rived from the 60 /xm luminosity, and 



SFRi.4GHz = 5.5 X 



i 1.4GHz 



4.6 X 102i(WHz-i) 



(2) 



(Haarsma et al. 2000; Condon 1992) for the SFR 
from the 1.4 GHz luminosity. In both cases the 
factor of 5.5 converts from the mass range 5 — 
100 Mq to the mass range 0.1 — 100 Mq, assuming 
a Salpeter IMF (Haarsma et al. 2000, for exam- 
ple). Also, Haarsma et al. (2000) point out that 
starburst galaxies may have IMFs that differ from 
the one adopted in Equation 2. In the worst case a 
maximum decrease of the SFR by a factor around 
5 is possible. Any changes to the slope of the IMF, 
however, should affect the radio and FIR luminosi- 
ties equally. The SFRs calculated for our BCD 
sample are given in Table 2 along with the dis- 
tances to the galaxies. The metallicities (oxygen 
abundance) also shown have been taken from IT. 
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HI masses have been derived using HI flux den- 
sities from the hterature with the distances given 
here, and are also shown in Table 2. 

Distance uncertainties (possibly up to a factor 
of a few) are the dominant source of measurement 
uncertainty for the SFRs. For this reason, the Fig- 
ures below do not include error bars on the SFRs. 
Histograms showing the distribution of the derived 
SFRs for our sample are presented in Figure 3. We 
also do not calculate SFR normalized to the op- 
tical area of the galaxy. This is firstly because of 
the complex optical morphologies for the objects 
in our sample (see Figure 2). Secondly, the galax- 
ies in our sample are essentially point sources for 
the 60 /xm and 1.4 GHz detections. Instead, we in- 
vestigate below the properties not only of the total 
SFR, but also of the "specific SFR" by which we 
mean the SFR normalized by HI mass. This is 
a distance-independent quantity, and a probe of 
relative star formation that should be reasonably 
independent of galaxy size or mass. Details of the 
specific SFRs and HI masses are given in Section 4 
below. 

The above calibrations of SFR from luminosity 
are derived independently of metallicity, although 
it is worth noting that the 1.4 GHz calibration is 
based on the observed Galactic relation between 
supernova rate and nonthermal radio luminosity. 
An implicit assumption in the case of estimating 
SFR from the 1.4 GHz luminosity is that there is 
no contribution from an active nucleus (Seyfert 
or LINER activity) to the radio emission. This 
assumption seems to be justified for the sample of 
BCDs we are investigating, since all of them have 
Hll-like spectra. 

4. Results 

The NVSS catalog yields flux densities for 
a significant fraction of our sample. Table 2 
and Figure 3 show that the SFRs derived from 
NVSS fluxes cover almost five orders of magni- 
tude. Clearly, there are real differences in BCD 
SFRs, since the range seen is much larger than 
the potential distance uncertainties. From the 
23 NVSS-detected objects we calculate a median 
SFR of 0.27MQyr-i. This compares weh with 
the median SFR, 0.36 Mgyr"!, of the 18 IRAS- 
dctcctcd BCDs. The IT sample of BCDs being 
examined, however, is likely closer to being a flux- 



limited than a volume-limited sample, favoring 
more luminous galaxies with higher SFRs, and 
the NVSS sensitivity limit strengthens this bias. 
Given that the radio and FIR detections arc sub- 
ject to such Malmquist bias, artificially raising 
the median value, the median SFR quoted above 
should not be over-interpreted. To provide some 
perspective, an approximate uncertainty for this 
median can be estimated by examining the range 
of SFRs spanned by the 3 (~ \/TT) galaxies either 
side of Mrk 1416 (having the median SFR= 0.27). 
This gives a range of an order of magnitude, from 
0.13 to 1.3MQyr-i. 

Figure 4 shows the comparison between SFRs 
derived from 1.4 GHz (FIRST and NVSS) and 
60 fxm luminosities. A good agreement exists be- 
tween the NVSS-based 1.4 GHz SFRs and the 
IRAS-based 60 /um SFRs. It is evident from Ta- 
ble 1 and this Figure that the FIRST flux den- 
sities tend to be systematically lower than those 
from the NVSS for most of the galaxies detected 
by both surveys. This is also true for many of 
the FIRST upper limits where an NVSS detec- 
tion exists. The probable reason for this is the 
difference between the two surveys in their sensi- 
tivity to extended emission. FIRST, made using 
the VLAs B-array, has much higher angular reso- 
lution than NVSS, made with the more compact 
D and DnC arrays, but at the expense of sensitiv- 
ity to more extended emission. FIRST images are 
typically only sensitive to structure extended over 
less than about 2', while the NVSS images will de- 
tect emission from sources up to about 15' across 
(Perley 2000). As a result the NVSS flux density 
measurements tend to be somewhat larger than 
those from FIRST, since many of the galaxies in 
our sample are larger than the extent over which 
FIRST is most sensitive. In subsequent analysis 
where 1.4GHz-derived SFRs are used we thus rely 
on SFRs derived from NVSS measurements, both 
for this reason and also since more of our sample 
galaxies have NVSS detections. 

Figure 5 presents SFRi.4ghz as a function 
of distance. The NVSS detection threshold is 
marked. As may be expected from the tight rela- 
tion in Figure 4, a diagram showing SFRgo^m ver- 
sus distance is very similar. An interesting point 
to note from this Figure is the absence of any mod- 
erate to high SFR objects closer than « lOMpc, 
that is, in the region where Hubble Space Tele- 



5 



scope (HST) observations (with WFPC2 or NIC- 
MOS instruments) are judged able to resolve in- 
dividual stars within the target galaxies. This is 
mainly a result of high SFR objects having lower 
space-densities than low SFR galaxies, combined 
with the small volume being probed at this dis- 
tance. The implication here is that HST-bascd 
CMD investigations of nearby galaxies are re- 
stricted to lower-SFR objects, in good agreement 
with the low SFRs derived from CMD modeling. 

The 1.4 GHz flux on which SFRi.4ghz is based 
is thought to measure the present-day SFR on a 
similar timescale as the lifetime of massive stars 
(Condon 1992). We now explore the correlation 
between various chemical abundances from IT and 
absolute and specific SFR. The chemical elements 
that we emphasize are associated with progeni- 
tors of different masses, and this translates into 
different enrichment timescales. Specifically, oxy- 
gen is produced by massive stars and returned 
quickly to the interstellar medium (ISM) by way 
of supernova (SN) type II explosions, nitrogen is 
thought to be produced in intermediate-mass stars 
and returned to the ISM in the asymptotic giant 
branch (AGB) wind phase (but see IT for a con- 
trary view), while iron is considered to reflect the 
ISM's pollution by low mass stars through SNIa 
events, and should lag the most behind the oxygen 
enrichment. 

In the top panel of Figure 6 a weak positive 
trend of 12-hlog(0/H) (from IT) with absolute 
SFRi.4GHz can be seen. The Spearman rank test 
shows that the null hypothesis (no correlation be- 
tween 12+log(0/H) and log(SFR)), using only 
the detections, has a probability of only about 
0.6%. An unweighted linear regression yields 
12-hlog(0/H) = 0.091og(SFR) + 7.90 (with the 
uncertainty of the slope being (Ja = 0.04 and the 
uncertainty of the intercept being ai, — 0.05). 
The extremely metal-poor galaxies I Zw 18 and 
SBS 0335—052 lie well away from the trend shown 
by the majority of the sample, effectively lower- 
ing the intercept. If we exclude I Zw 18 and 
SBS 0335—052 from the sample, than the Spear- 
man rank correlation coefficient increases slightly, 
and the probability of a chance correlation drops 
to just below 0.5%, while the hnear regression 
gives 12-hlog(0/H) = 0.091og(SFR) + 7.96 (with 
CTa = 0.03 and db = 0.03). It is this linear fit 
that is shown in Figure 6 as a solid line, with the 



dashed lines reflecting the error in the intercept. 
(Notice the fit results reported here are slightly 
different from the preliminary results of Schulte- 
Ladbeck et al. (2001b) because here wc use more 
up-to-date values for the Hubble constant and the 
Virgo distance.) We conclude that BCDs with 
high present-day total SFRs clearly also have high 
oxygen abundances. 

If we replace absolute SFRs with specific SFRs 
(i.e. SFR/Mhi) a weak positive trend seems to re- 
main (HI masses are discussed below). The Spear- 
man rank test, however, indicates the hypothesis 
that 12+log(0/H) and log(SFR/MHi) are corre- 
lated has a 15% probability of occur ing by chance. 
As a result it seems there is little correlation, if 
any, between oxygen abundance and SFR/Mhi- 

In the middle panel of Figure 6 a weak positive 
trend of log(N/0) (from IT) with absolute SFR 
can be seen. The extremely metal-poor galaxies 
I Zw 18 and SBS 0335-052 do not distinguish 
themselves on this plot, lying well within the dis- 
tribution of values. The Spearman rank test shows 
that the probability for a chance correlation is 
about 0.7%. A linear regression yields log(N/0) 
= 0.091og(SFR) - 1.42 (with aa = 0.02 and cxfe 
= 0.03). Conversely, the Spearman rank test for 
log(N/0) with the logarithm of the specific SFR, 
log(SFR/MHi), shows the probability for a chance 
correlation is 5%. 

Lastly, we investigate the dependence bewteen 
[0/Fe] (from IT) and log(SFR) or log(SFR/MHi) 
(see bottom panel of Figure 6). The Spearman 
rank test gives such small correlation coefficients 
that the data are likely to be uncorrected. The 
linear regression of [0/Fe] on log(SFR) results in 
a constant, [0/Fe] = 0.43 (with cja = 0.04 and C7(, 
= 0.04). Figure 6 shows that SBS 0335-052 is 
quite removed from the rest of the sample BCDs, 
but I Zw 18 lies well within the distribution. 

Figure 7 explores the dependence of various 
quantities on galaxy mass. The width of the HI 
emission line profile at 20% of peak, W20, is taken 
as an indicator of galaxy total mass. HI pro- 
file widths are often used to calculate galaxy to- 
tal masses, but require radii and inclinations. As 
can be seen from Figure 2, radii and inclinations 
are difficult to judge for BCDs because the star- 
forming regions dominate their appearance. W20, 
on the other hand, also may not be a solid indica- 
tor for rotational velocity if turbulent motions are 
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significant. 

With these caveats in mind, we plot in the up- 
per panel of Figure 7 the HI mass versus the W20 
velocity. HI and W20 data were taken from Thuan 
et al. (1999a) (circles) or Thuan & Martin (1981) 
(triangles). The Spearman rank correlation coeffi- 
cient of log(MHi) and log(W2o) is high, indicating 
a less than 0.05% probability that this is a chance 
correlation. In other words, we can be confident 
that HI mass increases as total mass increases. A 
linear regression gives log(MHi) = 3.61og(W2o) + 
1.3 (with (Ta = 0.6 and Ub = 1-3), and is shown in 
Figure 7. 

The Spearman rank test for total SFR and to- 
tal mass rejects the null hypothesis; correlation by 
chance has about a 0.5% probability (middle panel 
of Figure 7). The linear regression is log(SFR) 
= 4.81og(W2o) - 10.5 (with aa = 1.0 and db = 
2.0). The higher SFRs are clearly being hosted by 
higher mass galaxies, although there is a signifi- 
cant scatter in the trend (for a given rotational ve- 
locity or total mass, the SFR may vary over more 
than an order of magnitude). 

Finally, we also investigate the correlation be- 
tween metallicity and total mass. The lower panel 
of Figure 7 shows that the lowest metallicity ob- 
jects, such as I Zw 18 or SBS 0335—052, do not 
have the lowest rotational velocities. Indeed the 
Spearman rank test indicates that the probabil- 
ity of a chance correlation between 12-|-log(0/H) 
and log(W2o) is quite high, almost 10%. A higher 
total mass does not necessitate a higher oxygen 
abundance. 

Figure 8 shows the variation of SFRi.4ghz with 
HI mass. The data were taken from Thuan et 
al. (1999a) (circles) and from Thuan & Martin 
(1981) (triangles). Because of our interest in ex- 
tremely low-mctallicity objects, and because it is 
one of the BCDs with a TRGB distance, we added 
UGC 04483 to the Figure, using the HI flux from 
Huchtmeier & Richter (1986). This data point is 
shown as a square. The higher SFRs are clearly 
being hosted by galaxies with higher HI masses 
although there is a significant scatter in the trend 
(again, as with W20, for a given HI mass the 
SFR may vary over more than an order of magni- 
tude). This scatter is reflected in the wide range 
of timescales for gas consumption, from 0.3 Gyr to 
more than 10 Gyr, indicated by the solid lines in 
Figure 8. I Zw 18 and SBS 0335-052 do not stand 



out in any way from the rest of the sample in this 
picture. 

The results can be summarized as follows: 

1. The median SFR in the IT sample is 

0. 3M0yr~^ (with an uncertainty up to a 

factor of three or so, and it may be an upper 
limit due to Malmquist bias). 

2. Total SFRs in the IT sample range over at 
least four orders of magnitude. 

3. SFRs calculated from IRAS and from NVSS 
data are consistent with SFReo^m = SFRi.4ghz, 

1. e., these galaxies obey the FIR/radio cor- 
relation valid for large spirals. 

4. While BCDs in this sample are close, z<0.06, 
there is still a lack of high-SFR galaxies 
within the range of HST single-star photom- 
etry (about lOMpc). This is because such 
objects have space densities low enough that 
they are unlikely to appear in the small vol- 
ume probed out to this distance. 

5. We find a high probability of a positive corre- 
lation between 12 + log(0/H) or log(N/0) 
and total SFR, but [0/Fe] and total SFR 
appear to be uncorrclatcd. For all three 
element abundances, the probability of 
a correlation with SFR/Mhi (a distance- 
independent quantity) is low. 

6. Taking W20 as an indicator of total mass, 
we find a high probability for a positive cor- 
relation between HI mass and total mass, 
and between total SFR and total mass. The 
probability for a correlation of metallicity 
and total mass, on the other hand, is small. 

7. The "typical" timcscalc for gas consumption 
in the IT BCD sample is a few Gyrs. None of 
these galaxies is in danger of running out of 
gas in the very near future. However, these 
BCDs will die sooner than the typical dwarf 
Irregular galaxy. 

5. Discussion 

5.1. The Far-Infrared/Radio Correlation 
for BCDs 

We find excellent agreement between SFRgo^m 
and SFRi.4GHz using the calibration equations of 
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Condon (1992), Cram ct al. (1998), and Haarsma 
et al. (2000). Specifically, and as shown by Fig- 
ure 4, the data are consistent with SFReo/xm = 
SFRi.4GHz- This consistency, however, does not 
necessarily imply that the correct SFR has been 
determined. As an independent check of the SFRs 
from FIR and radio luminosities we compare them 
with results from modeling of color-magnitude di- 
agrams. CMD studies have their own but dif- 
ferent set of limitations, in terms of sensitiv- 
ity, color biasing, and fractional galaxy area ob- 
served. We have only two cases in the IT sample 
for which there arc both distances and present- 
day SFRs available from single-star photometry. 
For VII Zw 403, SFRgo^m = 0.0096 Moyr"! 
and SFRi.4GHz = 0.004 Moyr"!, while CMD 
modehng of HST/WFPC2 V and I photome- 
try yields O.Oll±O.OO8M0yr-i (Crone et al. 
2002). For Mrk 209 (I Zw 36) the SFRi.4ghz 
is 0.0217 M0yr-i, whereas HST/NICMOS J 
and H data were modeled to yield a SFR of 
0.025 ± 0.004 Mgyr-i (Schulte-Ladbcck et al. 
2001a). We conclude that CMD modeling vali- 
dates the FIR and radio SFRs to the extent that 
it is possible with the limited amount of data. 

Astrophysical interpretations of the FIR/radio 
correlation suppose that both types of emission 
have their origin in newly forming massive stars. 
The FIR emission originates from the heating of 
dust by UV radiation, while the radio emission 
is produced by synchrotron radiation from cos- 
mic ray electrons accelerated by shocks associ- 
ated with supernova remnants (Xu 1990). Both 
SFReo/jm and SFRi.4ghz are considered to be mea- 
sures of the present-day SFR. The synchrotron 
component of the radio emission is thought to 
dominate for a star-forming galaxy observed at 
1.4 GHz, with a much smaller contribution from 
the thermal Bremsstrahlung of the HII regions. 
With a typical synchrotron spectrum for the non- 
thermal component {Si, oc i^""'^), less than 1/8 of 
the total 1.4 GHz flux is thermal (Condon 1992). 
Equation 2 above makes the assumption that su- 
pernovae generate similar amounts of radio lumi- 
nosity to that observed in our own Galaxy, and 
the detection of any nonthermal radio luminosity 
in external star-forming galaxies implies the pres- 
ence of magnetic fields. 

The calibration above for SFRi.4ghz differs 
from the calibration given by Izotova et al. (2000). 



Izotova et al. studied SFRs in a sample of BCDs 
drawn from the Second Byurakan survey, but 
they find SFRi.4ghz ~ 3.5 SFReo/xm- Their cal- 
ibration for SFRi.4ghz from 1.4 GHz luminosity 
makes the assumption that the radio emission is 
purely thermal, although they acknowledge that 
a significant part of the 1.4 GHz radio emission 
should be nonthermal emission related to super- 
nova events. This explains their overestimation of 
1.4 GHz luminosity derived SFRs. The calibration 
for SFRgo/jm given by Izotova et al. also differs 
slightly from the one adopted above (but see dis- 
cussion of SFRfir calibrations given by Kcnnicutt 
1998). It can easily be shown that the objects in 
the sample of Izotova et al. follow the FIR/radio 
correlation, and indeed the SFRs derived for their 
sample by using the calibrations of Equations 1 
and 2 result in highly consistent SFRs from both 
60 iim and 1.4 GHz luminosities. 

There axe suggestions that the thermal compo- 
nent of the radio luminosity in BCDs is larger than 
that typical for normal spirals. Klein et al. (1991) 
found the radio spectra of BCDs are generally flat- 
ter than those of spiral galaxies. They comment 
that they find it surprising that BCDs do follow 
the same FIR/radio correlation as do normal spi- 
rals, given the strongly varying SFRs in BCDs 
(further discussed below) and the lack of cosmic- 
ray confinement. They suggest that a lower ratio 
of nonthermal-to-thermal radio emission in BCDs 
could be offset by a similar deficiency in FIR emis- 
sion. Deeg et al. (1993) re-observed a few of 
these BCDs and also found a variety of shapes for 
their radio spectra, but conclude that a nonther- 
mal component is present in most. Clearly, these 
data suggest that the radio spectra of BCDs at 
1.4 GHz have differing amounts of nonthermal-to- 
thermal emission, hence it is indeed surprising that 
our investigation yields identical FIR and 1.4 GHz 
SFRs when the calibrations of Equations 1 and 2 
are applied. There are no linear radio polarization 
observations of BCDs in the literature, and there- 
fore we have no direct measurements of their large- 
scale magnetic fields. The presence of a one-to-one 
FIR-radio correlation in our BCD sample suggests 
that magnetic fields in these systems need to be of 
at least sufficient strength to extract synchrotron 
radiation from relativistic electrons accelerated by 
supernova remnants. 

The conversion of 60/tim luminosity to SFRgo^m 
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rests on the assumption that dust absorbs well in 
the UV, and that it re-radiates effectively in the 
IR, so that the ultraviolet radiation of massive, 
young stars is absorbed and rc-cmittcd as FIR 
luminosity. Clearly the conversion factor must 
depend on the amount of dust present. There 
has been a persistent notion in the literature that 
metal-poor dwarf galaxies are virtually dust free. 
This is not the case. We emphasize that although 
our Table 2 shows no FIR detections for I Zw 18 
and SBS 0335—052, it is now known that both 
contain appreciable amounts of dust (Thuan et al. 
1999b; Dale ct al. 2001; Cannon et al. 2001). In 
any event, for the SFR60/tm/SFRi.4GHz ratio to be 
unity, it is not as much the total amount of dust in 
dwarfs that matters, but rather, that there is suf- 
ficient dust to produce a FIR luminosity that off- 
sets the radio luminosity. This idea supports the 
suggestion that the local properties of individual 
star-forming regions are more important in gen- 
erating the FIR/radio correlation than the global 
properties of the galaxies that host them. 

Is metallicity likely to affect the calibrations of 
SFR to luminosity? Lisenfeld & Ferrara (1998) 
discuss dust-to-gas ratios as a function of metal- 
licity in dwarf galaxies. Hirashita et al. (2001) 
incorporate metallicity into the conversion factor, 
and find it can vary by about a factor of 5 through 
differences in dust content (dust-to-gas ratios) be- 
tween our Galaxy and others. In order to inves- 
tigate whether there is any metallicity bias in the 
SFRs we derived for BCDs, we display in Fig- 
ure 9 SFRi.4GHz/SFR60;.m versus 12-hlog(0/H). 
This ratio spans a range of about 0.3 to 1.8 for 
metallicities of 12-hlog(0/H) = 7.7 to 8.3. As 
can be seen from this Figure we do not have mea- 
surements at both frequencies for a number of the 
extremely metal-poor objects. For three objects 
with a 60 /xm detection but a 1.4 GHz upper limit, 
we sec upper limits to this ratio between 0.8 and 
1.5. For six objects with 1.4 GHz detections but 
60 yum upper limits, this ratio shows lower lim- 
its ranging from 0.6 to 2.5. Both I Zw 18 and 
SBS 0335—052 fall into the latter category and 
have lower limits of 1.5. The mean metallicity of 
the detected sample is 12+log(0/H) = 7.98. IT 
use the scale on which the Sun has 12-|-log(0/H) = 
8.91, so we are effectively exploring the FIR/radio 
correlation at about 0.1 of the solar metallicity. 
In our sample, the mean SFRi.4GHz/SFR60/jm is 



1.02±0.40, identical to the expected value of 1. 
No trend with metallicity is seen. 

In starburst galaxies with highly episodic star- 
formation histories (SFH), a time-delay between 
the maximum of the FIR and the nonthermal ra- 
dio emission has been hypothesized (Lisenfeld et 
al. 1997). Thermal radiation and re-radiation re- 
spond instantaneously to the presence of newly 
born massive stars, whereas nonthermal radiation 
does not occur until the end of the life of massive 
stars, when they explode as SNe, some 10^ yrs 
later. If time-delays were important we would ex- 
pect a large scatter of SFRi.4GHz/SFR6o^m about 
the mean, depending on exactly when we observe 
a galaxy in its changing SFH. Indeed, BCDs might 
be suspected of this effect if they experience highly 
episodic starbursts, whereby very short bursts 
(of about 3.5 Myr, Mas-Hesse & Kunth 1999) 
are separated by very long quiescence periods of 
a Gyr or more (but see Schulte-Ladbeck et al. 
2001a). Weilbacher & Fritze-v.Alvensleben (2001) 
recently modeled the size of the offset in the Ha 
luminosity to SFRna calibration that would arise 
due to episodic starbursts. The conversion fac- 
tor ranges over two orders of magnitude for the 
benchmark SFHs that they explore. We observe 
a 1 CT variation of ±0.40 about the mean ratio for 
SFRi.4GHz/SFR6o^m, & much smaller scatter, al- 
though we admittedly are investigating an incom- 
plete sample. It is possible that time-delay effects 
contribute to this scatter about the mean, but they 
clearly do not dominate the conversion factor. 

To summarize, our finding of a metallicity- 
independent FIR/radio ratio of unity implies that 

the conversion factors that are currently being 
used to turn FIR and radio luminosity to SFR 
seem to remain valid at 0.1 of the solar metallicity. 
(In Paper II, we shall explore the conversion factor 
that translates Ha luminosity to SFRua for these 
objects.) We also infer that BCDs are required to 
have supernova rates, radio luminosities per SN, 
magnetic fields, and dust contents in some scaled 
proportion to those of the large spirals, to which 
the FIR and radio calibrations were originally tai- 
lored. 

5.2. SFR, Metallicity, and Mass 

A major difficulty for chemical evolution mod- 
els of BCDs has been to explain their low oxygen 
abundances in spite of their active star formation. 
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In this section we confront the data assembled here 
with current theoretical models for dwarf galaxy 
evolution. 

We see a correlation of HI gas mass with to- 
tal mass, and between total SFR and total mass, 
in the sense that galaxies with more total mass 
also have a higher gas mass, and support a higher 
present-day SFR. This comes as no surprise since 
the large range of SFRs among different galaxies 
is attributed in part to the large range; in galaxy 
masses (Kennicutt 1998). We see a broad trend 
of increasing SFR with increasing HI mass, but 
there is also a wide range of SFRs for a given HI 
mass. This is explained as being due to differ- 
ent gas-consumption timescales. Here we have a 
hint of some variation in the physical properties 
of BCDs. The gas depletion timescale of BCDs in 
the IT sample is typically a few Gyrs. This is an 
order of magnitude shorter than those in the dwarf 
Irregular sample studied by van Zee (2001). Our 
data also show a positive correlation between oxy- 
gen abundance and SFR, in the sense that galaxies 
with higher global SFRs also show higher oxygen 
abundance. This similarly comes as no surprise as 
the massive stars that are used to gauge the SFR 
are also the sites of oxygen nucleosynthesis. (How- 
ever, there is some as yet unsettled discussion in 
the literature as to whether the oxygen so mea- 
sured is related to the present generation of stars, 
or whether it stems from prior ones, e.g. Larsen et 
al. 2001 and references therein.) 

While SFR correlates with total mass and oxy- 
gen abundance correlates with SFR, there is no 
correlation between oxygen abundance and total 
niMss in our sample (see Figure 7) . We further find 
that SFR/Mhi versus HI mass shows little corre- 
lation, if any, and we also find no correlation be- 
tween oxygen abundance and SFR/Mhi- Lequeux 
et al. (1979) first presented a relationship between 
metallicity and total mass for a small sample of 
Irregular galaxies and two BCDs (I Zw 18 and 
II Zw 40). Skillman et al. (1989) later derived 
the well known metallicity-luminosity relation for 
Irregular galaxies. In Irregular galaxies, the oxy- 
gen abundance rises with absolute blue magni- 
tude of the host galaxy (but see Hildago-Gamez 
& Olofsson 1998; Hunter & Hoffmman 1999, for 
the opposing view). Note that Skillman et al. 
(1989) excluded BCDs from their sample. The 
metallicity-luminosity relation is frequently inter- 



preted as a metallicity-mass relation, to the extent 
that it is now being used to gauge galaxy assem- 
bly in the intermediate-redshift Universe (Kobul- 
nicky & Zaritsky 1999; CaroUo & Lilly 2001). In- 
deed, there is a general positive correlation across 
the Hubble sequence of oxygen abundance with 
Ms, that holds for over 10 orders of magnitude in 
Mb (Hunter & Hoffmman 1999). It is not as yet 
clear that a metallicity-luminosity relation holds 
for BCDs alone (Hunter & Hoffmman 1999, Fig- 
ure 7). Similarly, our limited data do not sup- 
port the existence of a mass-metallicity relation 
for BCDs. 

Since the seminal work of Dekel & Silk (1986), 
the extent to which mass loss affects the evolution 
of low-mass galaxies has attracted much attention. 
Such work is particularly interesting since low- 
mass systems are galactic building blocks in hier- 
archical scenarios of galaxy formation, and could 
contribute to the pollution of the intergalactic 
medium at early epochs. Recent models for dwarf 
galaxies have addressed the role of dark matter 
halos for the evolution of low-mass systems. Mac 
Low & Ferrara (1999) find that mass loss occurs in 
all galaxies with gas masses below w 10^ Mq, but 
complete loss of gas (blowaway) is possible only for 
galaxies with very low HI masses, less than a few 
times 10^ Mq. The HI masses in the IT sample 
range from several times 10^ Mq, to several times 
10^ M0. They match extremely well the HI mass 
range to which the models by Mac Low & Ferrara 
apply. This suggests that some mass loss occurs 
in all of the BCDs observed here. Mac Low & Fer- 
rara also point out that while the fraction of cool 
gas lost in dwarf galaxies is generally small, most 
of the metals mixed with hot gas are able to leave 
these galaxies. This helps to keep their metallici- 
ties low. Ferrara & Tolstoy (2000) re-iterate that 
even in the presence of a metal-enriched outflow, 
the predominant mechanism for gas consumption 
is the conversion of gas into stars, and not mass 
loss. 

Figure 10 shows the relation between oxygen 
abundance and gas (HI) mass for the IT sample of 
BCDs. Three dashed lines approximate the mod- 
els of Ferrara & Tolstoy (2000) for dark-to-visible 
(HI) mass ratios of 0, 10, and 300, respectively. 
Our Figure may be compared with the equiva- 
lent Figure for gas-rich Irregular galaxies in Fer- 
rara & Tolstoy (2000, their Figure 4). The data 
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for BCDs arc generally found to overlap with the 
model parameters. The BCD sample also shows 
considerable overlap with the "nearby dwarfs" and 
"low-surfacc-brightncss dwarfs" samples to which 
Ferrara & Tolstoy compared their models''. Our 
data thus lend additional support to the models of 
Mac Low & Ferrara and Ferrara & Tolstoy. This 
also means that most galaxies classed as BCDs are 
dark-matter dominated. 

Ferrara & Tolstoy predict total SFRs for a lim- 
ited choice of parameters. These SFRs are quite 
small, and allow a variation over a few orders of 
magnitude, owing to differences in dark-to-visible 
mass ratios and time evolution. The small global 
SFRs of BCDs that we measure here are broadly 
consistent with these predictions, although our 
median SFRs appear to be an order of magni- 
tude higher than their model predictions, as is 
the fact that we observe a wide range of SFRs. 
The effect of Malmquist bias on our sample may 
partially explain this result, biasing our median 
SFRs to higher values. Ferrara & Tolstoy empha- 
size that there is a considerable spread in oxygen 
abundances for a given HI mass, and that this can 
be understood because higher dark-to- visible mass 
ratios cause a stronger gas compression, increase 
the SFRs, and in doing so, also increase the oxygen 
abundances. A large difference in dark-to-visible 
mass ratios among the BCDs can help us to recon- 
cile the lack of correlation for oxygen abundance 
versus HI mass and oxygen abundance versus W20 
while, at the same time, SFR correlates with oxy- 
gen abundance and also shows a general trend with 
HI mass, albeit with much scatter. 

5.3. I Zw 18 and SBS 0335-052 

I Zw 18 and SBS 0335-052 once again stand 
out from the rest of the BCDs in Figure 10. They 
actually fall in a "forbidden" area of the model 
parameter space. The fact that the two well- 
studied extremely metal poor systems lie outside 
these models (i.e., in the unphysical region below a 
dark-to- visible mass ratio of zero) reinforces their 
unusual properties. This discrepancy may be ex- 
plainable through uncertainties in the HI masses, 
as, for I Zw 18 and SBS 0335-052, these only need 
to be lowered by 1-2 orders of magnitude to move 

'^A solar oxygen abundance of 12 + log{0/H) 8.9, consis- 
tent with that used by IT, is assumed for all three samples. 



them close to dark-matter dominated terrain in 
this Figure. This may be feasible because both 
I Zw 18 and SBS 0335-052 have complex opti- 
cal and radio morphologies (van Zee ct al. 1998; 
Lequeux & Viallefond 1980; Pustilnik et al. 2001). 
If the applicable HI mass is that of the "main 
body", then we can lower the HI flux by a fac- 
tor of 2 or so. A smaller distance will also lower 
the HI mass (although in the case of I Zw 18 it 
could not be lowered by more than a factor of 2-3, 
since otherwise it should be resolved more effec- 
tively into individual stars with HST). Indepen- 
dent data on dark-to-visiblc mass ratios arc also 
somewhat uncertain due to the many assumptions 
involved. While in the case of I Zw 18, radio ob- 
servations support a dark-to-visiblc mass ratio of 
about 10 (van Zee et al. 1998; Lequeux & Vialle- 
fond 1980), the dynamical mass of SBS 0335-052 
is only about 2 5 times that in gas and stars 
(Pustilnik et al. 2001). In any event, I Zw 18 
and SBS 0335—052, which are unlikely to be form- 
ing their first stars now (Schultc-Ladbcck ct al. 
2000) yet show oxygen abundances close to the 
"floor" of 12-|-log(0/H)«7.2 allowed by the mod- 
els of Ferrara & Tolstoy (2000), could be inter- 
preted as galaxies that have a maximum time- 
averaged mass-loss rate (since their masses are a 
few times 10^ M©) and are thus unable to retain 
metals effectively. 

Figure 4 of Ferrara & Tolstoy contains two 
galaxies, UGC 04483 and Leo A, with apparently 
similar properties to I Zw 18 in terms of oxygen 
abundance, HI mass, and dark-to-visible mass ra- 
tio. UGC 04483 is also present in the BCD sample 
of IT used here. While Ferrara & Tolstoy show it 
in the "forbidden" region, the smaller distance de- 
rived in the meantime from single-star photometry 
(Dolphin et al. 2001), has resulted in a downward 
revision of its HI mass. This, combined with the 
0.2 dex higher oxygen abundance of IT, has moved 
UGC 04483 well into the dark-matter dominated 
regime in our Figure 10. 

Ferrara & Tolstoy (2000) present a case study 
of the SFH and the oxygen enrichment of Leo A. 
It will be interesting to see whether their model 
also agrees with the broad parameters of I Zw 18 
(and SBS 0335-052). One clear difference be- 
tween Leo A and I Zw 18 (or SBS 0335-052), is 
that the present-day SFR of the former is much 
smaller, by about 3 (4) orders of magnitude, than 
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the latter. 

It is worth recalhng the location of I Zw 18 and 
SBS 0335-052 on Figure 6. They appear to be too 
metal-poor for their SFR, or too intensely star- 
forming for their oxygen abundance, compared to 
the rest of the sample. Explaining this offset as 
an error in the measured SFR is unlikely, requir- 
ing the actual SFR of I Zw 18 (for example) to 
be more than 7 orders of magnitude less than es- 
timated here. The metallicity of I Zw 18 has been 
the subject of many observations and extensive 
measurement, and explaining the offset as an error 
in the measured metallicity is similarly unlikely. 

It will be interesting to sec whether the models 
of Ferrara & Tolstoy can reproduce the high SFRs, 
high gas-mass fractions, and low chemical abun- 
dances of I Zw 18 and SBS 0335-052. Consid- 
ering the optical and HI morphologies of I Zw 18 
and SBS 0335—052, it is possible that these galax- 
ies are involved in an interaction. The models of 
Ferrara & Tolstoy do not apply to merger situa- 
tions. 

6. Conclusions 

We investigated the SFRs of a small (50 ob- 
jects) but well-defined sample of BCDs, that of 
Izotov & Thuan (1999). One advantage of study- 
ing this sample is that ionized-gas abundances are 
readily available from the work of this team. In 
addition, HI fluxes and profiles arc in hand for 
many of the BCDs from the work of Thuan and 
collaborators. This allows us to investigate cor- 
relations between SFR, ionized gas chemistry, HI 
mass, and total mass. The main results are sum- 
marized at the end of Section 4. Here, we draw 
some conclusions about some of the results that 
we think have a wider implication. 

We find that SFRgo^m agrees with SFRi.4ghz 
in BCDs. This agreement exists over a wide range 
of SFRs (nearly five orders of magnitude in the 
sample that we investigate) , at low metallicity (on 
average 0.1 of solar), and at small galaxy masses 
compared to those of normal spirals. Unless sev- 
eral factors have conspired to yield SFRgo^m = 
SFRi.4GHz in our sample, it is logical to conclude 
that the calibration relations and their underly- 
ing assumptions are valid for BCDs. Furthermore, 
the existence of this type of scaling relation for a 
wide range of galaxies with difi'erent masses and 



sizes (from spiral galaxies to late-type dwarfs), 
and with a wide range of interactions (from merg- 
ing to isolated), suggests that in terms of the 
FIR/radio correlation, the local properties of in- 
dividual star-forming regions are more important 
than the global properties of the galaxies that host 
them. 

While BCD studies have benefitted greatly 
from recent advances in far-infrared astronomy, 
and there is now no doubt that even very metal- 
poor galaxies living in the present epoch contain 
dust, very little is known as yet about their radio 
properties. Radio polarization observations are 
non-existent but are needed to elucidate the na- 
ture of magnetic fields in BCDs. Magnetic field 
orientations, therefore are not well known. Be- 
cause the FIR/radio correlation, valid for large 
spiral galaxies, extends to BCDs, the presence of 
magnetic fields at some level seems highly likely. 
The role of magnetic fields is still poorly under- 
stood, however, and may contribute to the ob- 
served scatter in the FIR/radio correlation. 

Color-magnitude diagrams of nearby BCDs can 
serve as an external check on the validity of the 
SFRs derived from FIR or radio luminosities. The 
relevant data are difficult to obtain, since they re- 
quire very high spatial resolution, and while the 
few existing data indicate encouraging agreement, 
additional CMD studies are required. Single-star 
photometry is, and in the near future will remain, 
biased toward nearby objects with small SFRs due 
to instrumental limitations. In Paper II, we shall 
therefore investigate the SFRs of BCDs in the IT 
sample using the Ha luminosities as SFR indica- 
tors. 

The small scatter in SFRs about the mean 
FIR/radio ratio supports results of CMD models 
that show that BCDs do not have highly episodic 
SFHs (Schulte-Ladbeck et al. 2001a). Our results 
also indicate that the median SFRs of BCDs are 
somewhat smaller than has been deemed typical 
in the past. The gas consumption timescales are 
at a minimum a few hundred Myr, and more typ- 
ically they are of the order of a few Gyrs. In this 
respect, BCDs appear closer to death than dwarf 
Irregulars (compare, for example, with the recent 
depletion time scales in van Zee 2001). 

We find that there is no mass-metallicity rela- 
tion among the BCDs in the IT sample. The ab- 
sence of a mass-metallicity relation for BCDs can 
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be understood with the help of differential winds 
(Marconi at al. 1994). The more recent models of 
Mac Low & Ferrara (1999) and Ferrara fc Tolstoy 
(2000) , which incorporate dark matter and are ap- 
plicable to the HI mass range that we observe in 
BCDs, appear to predict well the global properties 
of BCDs, at least to first order. An interesting 
prediction of these models is that all galaxies with 
HI masses below 10^ M© are prone to lose mass 
through the feedback of star formation. Interest- 
ingly, while there is not a lot of cool gas lost, the 
loss of metals becomes more and more important 
the lower the HI mass of the galaxy (with a peak 
in the metal-loss rate at an HI mass of a few times 
10^ Mq). This indicates a progressively stronger 
violation of the assumption of a closed-box model 
for chemical evolution toward low galaxy masses. 
We surmise that the absence of a mass-metallicity 
relation in BCDs reflects just this metal-enriched 
mass-loss. The basis for the mass-metallicity re- 
lation breaks down below a gas mass of about 
109 Mq. 

I Zw 18 and SBS 0335-052 continue to stand 
out as extreme examples of BCDs whenever their 
chemical properties are being investigated. In 
terms of total SFR, HI mass, gas-consumption 
timescale and W20, they fit in well with other 
BCDs. But in terms of oxygen abundance ver- 
sus SFR and HI mass, they are extremes com- 
pared to the rest of the sample. I Zw 18 and 
SBS 0335—052 have gas masses of a few times 
10^ Mq where dwarf galaxies, according to the 
models of Ferrara & Tolstoy, are particularly vul- 
nerable to the loss of heavy elements. A detailed 
prediction of the SFH and the chemical enrich- 
ment of I Zw 18 based on the work of Ferrara & 
Tolstoy woiild bo particularly interesting as a pre- 
cursor to the results of future deep single-star pho- 
tometry. Chemo-dynamical models with "gasp- 
ing," rather than bursting, SFHs for I Zw 18 are 
already in work by Recchi et al. (2002). ("Gasp- 
ing" in this sense means long episodes of moderate 
star-forming activity separated by short quiescent 
intervals, see e.g., Tosi et al. 2002.) It will be 
interesting to see what physical insights can be 
gained from their simulations. 



Acknowledgments 

AMH gratefully acknowledges support provided 
by NASA through Hubble Fellowship grant HST- 
HF-01140.01-A awarded by the Space Telescope 
Science Institute, which is operated by the Asso- 
ciation of Universities for Research in Astronomy, 
Inc., for NASA, under contract NAS 5-26555. 

REFERENCES 

Becker, R. H., White, R. L., Helfand, D. J. 1995, 
ApJ, 450, 559 

Cannon, J. M., Skillman, E. D., Garnett, D. R., 
Dufour, R. J. 2001, AAS, 198 #09.14 

Carollo, C. M., Lilly, S. J. 2001, ApJ, 584, L153 

Condon, J. J. 1992, ARA&A, 30, 575 

Condon, J. J., Cotton, W. D., Greisen, E. W., Yin, 
Q. F., Pcrley, R. A., Taylor, G. B., Broderick, 
J. J. 1998, AJ, 115, 1693 

Cram, L., Hopkins, A., Mobasher, B., Rowan- 
Robinson, M. 1998, ApJ, 507, 155 

Crone, M., Schulte-Ladbeck, R. E., Greggio, L., 
Hopp, U. 2002, ApJ, 567, 258 

Dale, D. A., Helou, G., Neugebauer, G., Soifer, B. 
T., Frayer, D. T., Condon, J. J. 2001, AJ, 122, 
1736 

Deeg, H.-J., Brinks, E., Duric, N., Klein, U., Skill- 
man, E. 1993, ApJ, 410, 626 

Dekel, A., Silk, J. 1986, ApJ, 303, 39 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, 
J. R., Buta, R. J., Paturel, G., Fouque, P. 1991, 
"Third reference catalogue of Bright galaxies," 
(New York: Springer- Verlag) 

Dolphin, A. E., Makarova, L., Karachentscv, I. 

D. , Karachentseva, V. E., Geisler, D., Grebel, 

E. K., Guhathakurta, P., Hodge, P. W., Sara- 
jedini, A., Seitzer, P. 2001, MNRAS, 324, 249 

Fanclh, M. N., O'Connell, R. W., Thuan, T. X. 
1988, ApJ, 334, 665 

Ferrara, A., Tolstoy, E. 2000, MNRAS,313, 291 



13 



Frccdman, W. L., Madorc, B. F., Gibson, B. K., 
Ferrarese, L., Kelson, D. D., Sakai, S., Mould, 
J. R., Kennicutt, R. C, Jr., Ford, H. C, Gra- 
ham, J. A.,Huclira, J. P., Hughes, S. M. G., 
Illingworth, G. D., Maori, L. M., Stetson, P. B. 
2001, ApJ, 553, 47 

Haarsma, D. B., Partridge, R. B., Windhorst, R. 
A., Richards, E. A. 2000, ApJ, 544, 641 

Hildago-Gamez, A. M., Olofsson K. 1998, A&A, 
334, 45 

Hirashita, H., Inoue, A. K., Kamaya, H., Shibai, 

H. 2001, A&A, 366, 83 

Hopkins A. M., Miller C. J., Connolly A. J., Gen- 
ovese C., Nichol R. C., Wasserman L. 2002, AJ, 
123, 1086 

Huchtmeier, W.K., Richter, O.G. 1986, A&AS, 63, 
323 

Hunter, D. A., Hoffmman, L. 1999, AJ, 117, 2789 

Izotova, I., Parnovskya, S., Izotov, Y. 2000, 

NewAR, 44, 283 

Izotov, Y. I., Thuan, T. X., Lipovetsky, V. A. 

1994, ApJ, 435, 647 

Izotov, Y. L, Thuan, T. X., Lipovetsky, V. A. 
1997, ApJS, 108, 1 

Izotov, Y. I., Thuan, T. X. 1998, ApJ, 500, 188 

Izotov, Y. I., Thuan, T. X. 1999, ApJ, 511, 639 
(IT) 

Kennicutt, R. C. 1998, ARA&A, 36, 189 

Klein, U., Weiland, H., Brinks, E. 1991, A&A 246, 
323 

Kobulnicky, H. A., Zaritsky, D. 1999, ApJ, 511, 
118 

Kraan-Korteweg, R. C. 1986, A&AS, 66, 255 

Kunth, D., Ostlin, G. 2000, A&AR, 10, 1 

Larsen, T. I., Sommer-Larsen, J., Pagel, B. E. J. 
2001, MNRAS, 323, 555 

Lequeux, J., Rayo, J. F., Serrano, A., Peimbert, 
M., Torres-Peimbert, S. 1979, A&A, 80, 155 



Lequeux, J., Viallefond, F. 1980, A&A, 91, 269 

Lisenfeld, U., Volk, H.J., Xu, C. 1997, 
RevMexAA, 6, 81 

Lisenfeld, U., Ferrara, A. 1998, ApJ, 496, 145 

Mac Low, M.-M., Ferrara, A. 1999, ApJ, 513, 142 

Marconi, G., Matteucci, F., Tosi, M. 1994, MN- 
RAS, 270, 35 

Mas-Hesse, M., Kunth, D. 1999, A&A, 349, 765 

Pcrlcy, R. A. 2000, "The VLA 
Observational Status Summary," 

(http://www.aoc.nrao.edu/vla/obstatus/vlas/vlas.html) 

Popescu, C. C., Hopp, U., Rosa, M. R. 1999, 
A&AS, 350, 414 

Pustilnik, S. A., Brinks, E., Thuan, T. X., 
Lipovetsky, V. A., Izotov, Y. 1. 2001, AJ, 121, 

1413 

Recchi, S., Matteucci, F., D'Ercole, A., Tosi, M. 
2002, A&A, in press (astro-ph/0112450) 

Sage, L. J., Salzer, J. J., Loose, H.-H., Henkel, C. 
1992, A&A, 265, 19 

Sargent, W. L. W., Searle, L. 1970, ApJ, 162, L155 

Schulte-Ladbeck, R. E., Hopp, U., Crone, M. M., 
Greggio, L. 1999, ApJ, 525, 709 

Schultc-Ladbcck, R. E., Hopp, U., Crone, M. M., 
Greggio, L. 2000, in "The First Stars" , eds. A. 
Weiss, T. Abel, V. Hill, (Springer) p. 293 

Schulte-Ladbeck, R. E., Hopp, U., Greggio, L., 
Crone, M. M., Drozdovsky, I.O. 2001a, AJ, 121, 
3007 

Schulte-Ladbeck, R.E., Hopp, U., Greggio, L., 
Crone, M.M., Drozdovsky, I.O., Hopkins, A.M. 
2001b, in Proc. of the 40*^* meeting of the Bonn- 
Bochum Graduicrtcn KoUcg "Dwarf Galaxies 
and their Environment", eds. K. S. de Boer, 
R.-J. Dettmar, U. Klein, Shaker Verlag, p. 83 

Skillman, E. D., Kennicutt, R. C, Hodge, P. W. 
1989, ApJ, 347, 875 

Thuan, T. X., Martin, G. E. 1981, ApJ, 247, 823 



14 



Thuan, T. X., Lipovctsky, V. A., Martin, J.-M., 
Pustilnik, S. A. 1999a, A&AS, 139, 1 

Thuan, T. X., Sauvage, M., Madden, S. 1999b, 
ApJ, 516, 783 

Tonry, J. L., Dressier, A., Blakeslee, J. P., Ajhar, 
E. A., Fletcher, A. B., Luppino, G. A., Metzger, 
M. R., Moore, C. B. 2001, ApJ, 546, 681 

Tosi, M., Grcggio, L., Annibah, F., Aloisi, A. 2002, 

astro-ph/0112268 

van Zee, L. 2001, AJ, 121, 2003 

van Zee, L., Westpfahl, D., Haynes, M. P. 1998, 
AJ, 115, 1000 

Weilbacher, P.M., Fritze-v.Alvensleben, U. 2001, 
A&AL, 373, 9 

Xu, C. 1990, ApJ, 365, L47 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.0. 



15 



Table 1 
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Table 1 — Continued 
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Table 2 

Distances. SFRs and metallicities for the local BCD galaxies. 
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Table 2 — Continued 



Name 


D 


iSi^iieO/im 




12+log(0/H) 


Mhi 




(Mpc) 


(Moyr-i) 


(Moyr-i) 
NVSS (FIRST) 




(10«Mq) 


SBS 1319+579 


36.7 




<0.56 


8.10 + 0.01 




SBS 1331+493 


15.1 




<0.078 (<0.030) 


7.82 + 0.05 


3.5 + 0.2 


MRK 1486 


150.6 




6.2 


7.88 + 0.01 




CG 0389 


15.3 




0.14 (<0.032) 


7.59 + 0.01 


2.7 + 0.2 


CG 0413 


94.5 




<3.0 (<1.2) 


7.75 + 0.01 




MRK 0475 


13.6 




<0.072 (<0.027) 


7.93 + 0.01 


0.066 + 0.02 


CG 1258 


199.4 




<12.5 (<5.6) 


7.99 + 0.06 




SBS 1533+574 


53.2 


0.94 


1.7 


8.01 + 0.12 




MRK 0487 


15.3 


0.092 


<0.092 


8.06 + 0.04 


0.80 + 0.14 


MRK 0930 


75.3 


9.1 


10.6 


8.06 + 0.03 


30.2 + 0.3 



19 




15 



5000 lO'* 1.5x10-* 

Vq [km/sec] 
I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



o 

c 10 h 
(U 5 - 

El. 







I I 



Lr- 



1 I ■ I-M I I 1-1 



50 100 150 200 250 

d [Mpc] 



Fig. 1. — Histograms showing the hehocentric velocities and Virgo-corrected distances for the sample. 
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Fig. 2. — Optical images from the Digitized Sky Survey (greyscale) with 1.4 GHz radio contours in white 

from the NVSS ovcrlaycd. The lowest contour level is 1 mJy in all cases apart from UGC 04483 (0.5 mJy) and 
MRK 1416 (0.8 mJy). For the 6 objects detected at 1.4 GHz by FIRST, contours are shown in black, with the 
lowest contour level being 0.5 mJy. Each image shown is about 4' x 4'. The 24 galaxies with NVSS or FIRST 
detections are shown in this Figure (continued overleaf). From left to right, top to bottom the galaxies shown 
above are: UM 311; UM 420; SBS 0335-052; MRK 1089; MRK 0071; MRK 1409; UGC 04483; MRK 1416; 
SBS 0926+606. 
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Fig. 2.- Continued. From left to right, top to bottom the galaxies shown above are: I Zw 18; CG 0798; 
MRK 0036; MRK 0162; VII Zw 403; UM 448; MRK 0750; UM 462; SBS 1222+614. 
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Fig. 2. Continued. From left to right, top to bottom the galaxies shown above are: MRK 0209; MRK 0059 
MRK 1486; CG 0389; SBS 1533+574; MRK 0930. 
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Fig. 3. — Histograms showing the range of derived SFRs for the sample. For the NVSS and FIRST (1.4 GHz) 
estimates, the inclusion of the upper limits as if they had been detections are indicated by the dashed 
histograms. This figure primarily indicates the broad range of SFRs exhibited in BCD galaxies. 
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Fig. 4.— Comparison of SFReo^ixm with SFRi.4ghz from both NVSS (circles) and FIRST (triangles). The 
dashed line is the relation SFRgo^jm = SFRi.4ghz- Note that the measurements of SFRi.4ghz derived from 
FIRST are, in most cases, systematically lower than those from NVSS. This is likely to be a result of NVSS 
having greater sensitivity to the same extended structure as that detected by IRAS. FIRST, although having 
much higher resolution than NVSS, detects notably less of this more extended emission. 
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Fig. 5. — SFRi.4GHz as a function of distance. The dashed line corresponds to the sensitivity hmit of the 
NVSS catalogue (equivalent to a 2niJy detection). The hashed area indicates the region (less than about 
10 Mpc) within which it is feasible to perform color-magnitude diagram analyses using HST observations. 
Note the absence of any local BCDs with SFRs greater than about 0.1 yr~^. 
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Fig. 6. — Emission-line based abundance indicators for BCDs from IT such as (12 + fog(0/H)), log (N/0) 
and [0/Fe] as a function of SFRi.4ghz- A least squares fit to the objects with measurements (not upper 
limits) of both SFR and abundance is shown, along with the 1 a uncertainties for the intercept of the fit. 
Uncertainties on the slope are omitted for clarity. 
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Fig. 7. — SFRi.4GHz as a function of W20, a rotational velocity estimator based on HI emission line width. 
The circles refer to data from (Thuan et al. 1999a) while the triangles us the data from (Thuan & Martin 
1981). The least squares fits to the detections and la uncertainties for the intercept are again shown. 
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Fig. 8. — SFRi.4GHz as a function of Mhi (Thuan et al. 1999a; Thuan & Martin 1981). The square indicates 
data for the metal-poor galaxy UGC 04483, using the HI flux and flux error of Huchtmeier & Richter (1986). 
Naive estimates of gas consumption time, calculated simply by dividing the gas mass by the SFR, are shown. 
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Fig. 9. — To investigate any possible bias between the 1.4 GHz and 60 /um SFR indicators, we explore the 
relationship between SFRi.4GHz/SFR60;:im and metallicity, 12 + log(0/H). The lines show the mean values 
(solid) and the ±1 a values (dashed) for each of the quantities. There does not seem to be any particular 
trend visible, supporting the use and self-consistency of the two SFR indicators. 
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Fig. 10. — Mctallicity, 12 + log(0/H), as a function of Mhi- Dashed lines show dark-to-visible mass ratios 
after Ferrara & Tolstoy (2000). Note the extremely metal poor galaxies I Zw 18 and SBS 0335—052 again 
lying well below the rest of the population. 
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